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Abstract 

We present a homodyne detection scheme to verify Bell's inequality on cor- 
related optical beams at the output of a nondegenerate parametric ampli- 
fier. Our approach is based on tomographic measurement of the joint de- 
tection probabilities, which allows high quantum efficiency at detectors. A 
self-homodyne scheme is suggested to simplify the experimental set-up. 
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I. INTRODUCTION 



In 1935 Einstein, Podolsky and Rosen []I| proved the incompatibility among three hy- 
potheses: 1) quantum mechanics is correct; 2) quantum mechanics is complete; 3) the 
following criterion of local reality holds: "If, without in any way disturbing a system, we 
can predict with certainty [...] the value of a physical quantity, then there exists an element 
of physical reality corresponding to this quantity." The paper opened a long and as yet 
unsettled debate about which one of the three hypotheses should be discarded. While Ein- 
stein suggested to abandon the completeness of quantum mechanics, Bohr refused the 
criterion of reality. The most important step forward in this debate was Bell's theorem of 
1965 0, which proved that there is an intrinsic incompatibility between the assumptions 1) 
and 3), namely the correctness of quantum mechanics and Einstein's "criterion of reality". 
In Bell's approach, a source produces a pair of correlated particles, which travel along oppo- 
site directions and impinge into two detectors. The two detectors measure two dichotomic 
observables A{a) and B{i3) respectively, a and (3 denoting experimental parameters which 
can be varied over different trials, typically the polarization/spin angle of detection at each 
apparatus. Assuming that each measurement outcome is determined by the experimental 
parameters a and f3 and by an "element of reality" or "hidden variable" A, Bell proved an 
inequality which holds for any theory that satisfies Einstein's "criterion of reality" , while it 
is violated by quantum mechanics. Such a fundamental inequality, which allows an experi- 
mental discrimination between local hidden-variable theories and quantum mechanics, has 
been the focus of interest in a number of experimental works [Q. 

Unfortunately, Bell's proof is based on two conditions which are difficult to achieve 
experimentally. The first is the feasibility of an experimental configuration yielding perfect 
correlation; the second is the possibility of approaching an ideal measurement, which itself 
does not add randomness to the outcome. Since 1969, attention was focused on improving 
the correlation of the source on one hand, and, on the other, on deriving more general 
inequalities that take into account detection quantum efficiency or circumvent the problem, 
however, at the cost of introducing supplementary hypotheses (see Refs. [^), as the well 
known "fair sampling" assumption. Anyhow it was clear also to the authors of the same 
Refs. that these assumptions are questionable, and, as a matter of fact, it was proved 
that in all performed experimental checks the results can be reproduced in the context of 
Einstein's assumptions if quantum efficiency of detectors is less than 82.3%. However, no 
experiment has yet succeeded in realizing such a high value of quantum efficiency. 

In a typical experiment the source emits a pair of correlated photons and two detectors 
separately check the presence of the two photons after polarizing filters at angles a and (3, 
respectively. Alternatively, one can use four photodetectors with polarizing beam splitters 
in front, with the advantage of checking through coincidence counts that photons come in 
pairs. Let us denote by Pa,f3 the joint probability of finding one photon at each detector with 
polarization angle a and f3, respectively. In terms of the correlation function 

C{a,(3) =Pa,(3+Pa,^-Pa,p-Pa,0 , (l) 

Bell's inequality writes as follows 

B{a,/3,a\/3') = \Cia,/3)-Cia,/3')\ + \Cia',f3') + Cia',/3)\ <2, (2) 
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a and /3 being the polarization angles orthogonal to a and ^ respectively. In this letter we 
propose a new kind of test for Bell's inequality based on homodyne tomography (for a 
review see Ref. 0). In our set-up the photodetectors are replaced by homodyne detectors, 
which along with the tomographic technique can be regarded as a black box for measuring 
the joint probabilities Po,/?. The main advantage of the tomographic test is that it allows 
using linear photodiodes with quantum efficiency ri higher than 90% |jTO|. On the other hand, 
the method works effectively even with t] as low as 70%, without the need of a "fair sampling" 
assumption, since all data are collected in a single experimental run. With respect to the 
customary homodyne technique, which in the present case would need many beam splitters 
and local oscillators (LO) that are coherent each other, the set-up is greatly simplified by 
using the recent self-homodyne technique [|lTl. Another advantage of self-homo dyning is the 
more efficient signal-LO mode-matching, with improved overall quantum efficiency. 



II. THE EXPERIMENTAL SET-UP 

The apparatus for generating the correlated beams is a x*^^'' nonlinear crystal, cut for 
Type-II phase-matching, acting as a nondegenerate optical parametric amplifier (NOPA). 
The NOPA is injected with excited coherent states (see Fig. 0) in modes c^, cj, c?^, (i| all 
with equal intensities and at the same frequency tuo, c and d denoting mode operators for 
the two different wave-vector directions, and | and representing vertical and horizontal 
polarization, respectively. The NOPA is pumped at the second harmonic loj^. At the output 
of the amplifier four photodetectors separately measure the intensities /aj, /ft^, Jq^, Jft^ of 
the mutual orthogonal polarization components of the fields propagating at different wave 
vectors. A narrow band of the photocurrent is selected, centered around frequency <^ cuq 
(typically ti^o is optical/infrared, whereas is a radio frequency). In the process of direct 
detection, the central modes C|^^ and d-^^^ beat with c^oi^ sidebands, thus playing the role 
of the LO of homodyne detectors. The four photocurrents la-^ , , la^ , yield the value 



of the quadratures of the four modes 



= --^(a^(+) + a^(-)) , s = {a,b}, vr = {^, J} , (3) 

where a^(±) and &7r(±) denote the sideband modes at frequency ujq ± Q, which are in 
the vacuum state at the input of the NOPA. The quadrature is defined by the operator 

= |(ae~*'^ + a^e*"^), where cj) is the relative phase between the signal and the LO. The 
value of the quadratures is used as input data for the tomographic measurement of the 
correlation function C{a,f3). The direction of polarizers {a, (3) in the experimental set-up 
does not need to be varied over different trials, because, as we will show in the following, 
such direction can be changed tomographically. 

We will now enter into details on the state at the output of the NOPA and on the 
tomographic detection. In terms of the field modes in Eq. (|^) the spontaneous down- 
conversion at the NOPA is described by the unitary evolution operator 

f/(^) = exp (a{6L + - h. c] , (4) 

where ^ = x^^^lL/c is a rescaled interaction time written in terms of the nonlinear suscep- 
tibility x*^^) of the medium, the crystal length L, the pump amplitude 7 and the speed c of 
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light in the medium, whereas represents a tunable phase that can be varied by rotating 
the crystal around the optical axis |12]. The state at the output of the NOPA writes as 
follows 



1^) = (1 - \A\') Y: E A"+"e^^"l^, n, m, m) ^ l^i,^) ® 1^3,4) , (5) 

n=0 m=0 

where A = ^/|^| tanh|^| and \i,l,m,n) represents the common eigenvector of the number 
operators of modes a|, b^, a^, with eigenvalues i, I, m and n, respectively. The average 
photon number per mode is given by iV = |Ap/(l — |Ap). The four-mode state vector in 
Eq. factorizes into a couple of twin beams IV'1,2) and |V^3,4), each one entangling a couple 
of spatially divergent modes (aj, and a^, respectively). 

Notice that conventional experiments, concerning a two-photon polarization-entangled 
state generated by spontaneous down-conversion, consider a four-mode entangled state which 
corresponds to keeping only the first-order terms of the sums in Eq. (^, and to ignoring 
the vacuum component, as only intensity correlations are usually measured. Here, on the 
contrary, we measure the joint probabilities on the state @ to test Bell's inequality through 
homodyne tomography, which yields the value of B{a, f3,a', (3') in Eq. (Q). 



III. TOMOGRAPHIC TEST OF BELL'S INEQUALITY 



The tomographic technique is a kind of universal detector, which can measure any ob- 
servable O of the field, by averaging a suitable "pattern" function 7l[0]{x, (p) over homodyne 
data X at random phase (p. The "pattern" function is obtained through the trace-rule 

n[d]{x, 0) = Tr lOK^ix - x^) 



(6) 



where K^{x) is a distribution given in Ref. [T^. For factorized many- mode operators O = 
Oi ® O2 ® ... ® On the pattern function is just the product of those corresponding to each 
single-mode operator Oi,...,On labelled by variables (xi, 0i), By linearity the 
pattern function is extended to generic many-mode operators. 

Now we consider which observables are involved in testing Bell's inequality @. Let us 
denote by qa,/3{h l^rn^n) the probability of having i,l,m,n photons in modes ai,b^,a^, 6j 
for the "rotated" state 



Ui,3{a) and U2,4{P) being the unitary operators 

Ui^3{a) = exp a (aja^ — '^I'^L) 
U24(3) = exp [(3 - ^jfot 



(7) 

(8) 
(9) 



The probabilities in Eq. (||) can be written as pa,/3 = Pa,/3(1,1), Pa,f3 = Pa,f3{0,0), pa,p 
p„,/3(0, 1), and Pa^p = Pa,/3(1, 0), with 



PaA^.m) 



qa,f3{n, 1 — m,l — n,m) 

^ POO) 



(10) 



4 



and {n,m = 0, 1}. The denominator P(l, 1) in Eq. (1^) represents the absolute probabihty 
of having at the output of the NOPA one photon in modes a|, and one photon in modes 
independently on the polarization, namely 

P{^A)= J2 g„,/3(n, 1 -m, 1 -n,m) . (11) 

n=0,l m=0,l 

Notice that our procedure does not need a fair sampling assumption, since we measure in 
only one run, both the numerator and the denominator of Eq. (ITUp, namely we do not have 
to collect auxiliary data to normalize probabilities. On the other hand, since we can exploit 
quantum efficiencies as high as rj = 90% or more, and the tomographic pattern functions 
already take into account r], we do not need supplementary hypothesis for it. 

The observables that correspond to probabilities qa,f3{h l,m,n) in Eqs. ([ToD and ( |TI1) are 
the following 

I i, I, m, n)a^p i, /, m, n\ = IJl^ia) IJl^iP) \ i, I, m, n) { i, /, m, n\ U2,i{(3) f/i,3(a) • (12) 

After a straightforward calculation using Eqs. (|TD|), (|TT|) and (0), one obtains that P(l, 1) 
is measured through the following average AV of homodyne data 

P(l, 1) = { {Kl Kl + Kl Kl) {Kl + Kl K^) } , (13) 

where denotes the diagonal (n = 0, 1) tomographic kernel function for mode j, namely 

Ki^{n\K,{x-x^^)\n) . (14) 

The probabilities in the numerator of Eq. (|TOD are given by the average of a lengthy ex- 
pression, which depends on both the diagonal terms (14) and the following off-diagonal 
terms 

Ki = {0\K,{x - x^Jll) , KL = {l\K,{x - x^JlO) = {Ki)* . (15) 

Here we report the final expression for C(a,/3) of Eq. (0) 

C{a, (3) = -— ^^v{[cos(2a) [kI - Kf) + sin(2a) [kI + K\ Kl) ] x 

[cos(2/?) (kI Kt - Kl K^) + sin(2/3) (kI Kt + KIkX)]}. (16) 

Caution must be taken in the estimation of the statistical error, because C{a,(3) — 
and thus B{a, (3,a', P') in Eq. (0) — are non linear averages (they are ratios of averages). 
The error is obtained from the variance calculated upon dividing the set of data into large 
statistical blocks. However, since the nonlinearity of B introduces a systematical error which 
is vanishingly small for increasingly larger sets of data, the estimated mean value of B is 
obtained from the full set of data, instead of averaging the mean value of blocks. 
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IV. NUMERICAL RESULTS 



We now present some numerical results obtained from Monte-Carlo simulations of the 
proposed experiment. For the simulation we use the theoretical homodyne probability per- 
taining to the state in Eq. (|) which, for each factor is given by 



2 exp 



TT 



+4A2)(rf2 +4A2) 



(17) 



where Xi {i = 1, 2, 3, 4) is the outcome of the homodyne measurement for quadrature of the 
i-th mode at phase (pi, and 



^(^»+^^)A, 4.„ = ^^^, K=^-r^- (18) 



At] 



In Fig. H we present the simulation results for B in Eq. (H) vs the phase ip in the state of 
Eq. d^). The full line represents the value of B in Eq. with the quantum theoretical 
value C{a,P) given by 

C{a,P) = cos V9 sin 2q; sin 2/3 — cos 2q; cos 2/3 . (19) 

Quantum efficiency rj = 85% has been used, nonetheless notice that for (p = n (corresponding 
to a maximum violation with respect to the classical bound 2), the obtained value is over 
10 a distant from the bound. By increasing the number of homodyne data, it is possible to 
obtain good results also for lower quantum efficiency. In fact, by increasing the number of 
data to 8-10^, a value of 5(0, fyr, f , |) = 2.834±0.268 has been obtained for A = 0.5, ip = tt, 
and 7] as low as 65%. This result is to be compared with the quantum theoretical value of 
2a/2. In Fig. ^ the results of the measurement of B, for different simulated experiments 
using the same number of data, are presented for different detector efficiencies t]. Notice 
how the error bars decrease versus rj. 

For an order of magnitude of the data acquisition rate in a real experiment, one can 
consider that in a typical self-homodyne set-up with a NOPA pumped by a 2"^^ harmonic of 
a Q-switched mode-locked Nd:YAG laser, the Q-switch and the mode-locker repetition rates 
are 10 kHz and 100 MHz, respectively. Typical time of the boxcar integration is 10 ns, so 
that one sample per pulse can be collected. In summary, 10^ data samples can be obtained 
in 10^ s. In such an experimental arrangement, for a Q-switch envelope of 200 ns, the shot 
noise can be reached by the peak amplitude of the 5-MHz low-pass-filtered photocurrent. 
For more detailed experimental parameters, the reader is referred to Ref. flf 



V. CONCLUSIONS 



In conclusion we have proposed a test of Bell's inequality, based on self-homodyne to- 
mography. The rather simple experimental apparatus is mainly composed of a NOPA crystal 
and four photodiodes. The experimental data are collected through a self-homodyne scheme 
and processed by the tomographic technique. No supplementary hypotheses are introduced. 
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a quantum efficiency rj as high as 90% is currently available, and, anyway, r] as low as 70% is 
tolerated for 10^-10^ experimental data. We have presented some numerical results based on 
Monte-Carlo simulations that confirm the feasibility of the experiment, showing violations 
of Bell's inequality for over 10 a with detector quantum efficiency 77 = 85%. 
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FIGURES 
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FIG. 1. Experimental set-up for the tomographic test of Bell's inequality. PBS and BS denote 
'polarizing beam splitter' and 'conventional beam splitter' respectively. Input radiation modes a|, 
b^, and 6| are in the vacuum state, while modes C|, c^, dp (at laser frequency loq) are 
in a coherent state. At the output of the nondegenerate parametric amplifier (NOPA) the four 
photocurrents I at radiofrequency O are measured, yielding the value of quadratures of the field 
modes a ^ , , and . The outcome quadratures are then used to reconstruct the probabilities 
of interest through quantum tomography. 
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FIG. 2. Plot of B{a, /5, a', (3') vs the phase ip in the state of Eq. @ for a simulated experiment. 
The shaded area represents the classical region for B. The parameters of the simulation are: 
a = 0; (5 = Ivr; a' = ^; = |; quantum efficiency r/ = 85%; average photon number per mode 
N = 0.5. A total number of 10^ homodyne data (divided into 20 statistical blocks) has been used. 
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FIG. 3. Plot of B{a, l3,a' , l3') vs the quantum efficiency of the detectors for a series of simu- 
lated experiments. The shaded area represents the classical region for B. The parameters of the 
simulations are: a = 0; f3 = |7r; a' = j; = ^; (/9 = vr; = 0.5. A total number of 6 • lO''' 
homodyne data (in 20 statistical blocks) has been used for each simulation. 
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